Quantum kinetic equation and universal conductance fluctuations in graphene 
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We analyze universal conductance fluctuations (UCF) in graphene in the framework of diagram- 
matic perturbation theory in the metallic regime. It is shown that strong inter-valley scattering 
lifts the valley degeneracy of electronic states, whereas at weak inter-valley scattering two valleys 
contribute independently such that the variance of UCF would be expected to show sample- and 
geometry-dependent behavior. 
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The unusual chiral properties of charge carriers in 
graphene ^, 2s, 5, 4] have recently received a lot of atten- 
tion. Several theories have been developed 0, H, 0i B 
interpreting observation of quantum interference effects 
in gra phe ne, such as weak localization magnetoresis- 
tance [ifl, [ill and the Josephson proximity effect in 
superconductor-graphene-superconductor junctions 1131 . 
Low-temperature magnetoresistance measurements [lOl . 
[isl . have shown universal conductance fluctuations 
(UCF) which appear to be robust over a wide range of 
electron concentrations and magnetic fields, and numeri- 
cal simulations of transport in monolayer graphene with 
charged disorder showed samplc-to-sample variation of 
conductance [iSj . In this Communication we analyze 
UCF using the same framework as the earlier weak lo- 
calization studies 0]- Specifically, we study UCF in 
graphene with various types of disorder in the fully de- 
veloped metallic regime (kpl ^ 1), using quantum ki- 
netic equation for diffusive transport in graphene and 
technique of semiclassical Keldysh functions. 

The transport in graphene is determined by the low- 
energy properties of charge carriers in the vicinity of cor- 
ners (if-points) of hexagonal Brillouin zone JJJ, called 
valleys. In the case of monolayer graphene, this can be 
described using the Hamiltonian [ITI . Il8l | , 



= t)Sp + /iw(p) + U(r), 



(1) 



Here, the basis of bi-spinors =[(/'K+,yi, 4>'k+,b, 0k_ b, 
0k_,a] characterizes electronic amplitudes on two crys- 
talline sublattices of graphene [A and B). Es and A;, 
{I = X, y, z) are 4x4 matrices in the valley and sublat- 
tice spaces [l^, introduced in Ref.j^. The momentum 
p = p(cos if, sin if) is defined with respect to the K-points 
[Iq . The first, 'Dirac' term in ([T]) determines an almost 
linear spectrum e = ±vp of electrons. The trigonal warp- 
ing term, /i^, takes into account a slight trigonal asym- 
metry of the Fermi line of graphene in one valley (such 
that e(K, p) ^ e(K, — p)), which will be treated below as 
a weak perturbation. Due to the time-reversal symmetry 
of the system the trigonal warping has opposite sign in 
K and K' valleys, e(K,p) = e(K',— p), which is taken 
into account by the valley matrix structure of H. The 



time-reversal-symmetric disorder, 

y(r) = iu(r) + ^E,AjM,j(r), 

consists of the potential Iu{r) due to remote Coulomb 
charges in, or on, the surface of the substrate (the unit 
matrix in the valley and sublattice space) and a generic 
term, which takes into account all possible symmetry- 
breaking local perturbations. The disorder is character- 
ized by {usi{r)us'i'{r')) = 5{r - r')6ss'Su'Wsi, and this 
determines the corresponding scattering rates r^^^ = 
Ss,s'Si,i'TrjWs,i/h (where 7 = kp / i'^TTvh) is the Fermi den- 
sity of states). After averaging over impurity configura- 
tions the scattering rates should preserve rotational sym- 
metry of graphene which means that t^i = Tyi = t±i and 
r,„ = Ts±. Anticipating a little, two scattering 
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rates, t^^ — 4t7 



2t~ and = 4 t7 
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' ±z ' ^'zz 'i — ^ '_L_L 

scribe valley block-diagonal and valley block-off-diagonal 
parts of the symmetry-breaking disorder potential re- 
spectively, whereas the total scattering rate is defined 
as T-i = Tq-i + T-i + 2tZI + 2t;1 + At^I. 

To characterize the UCF we evaluate the variance of 
conductance, {5Q'^) = {Q^) — {G)^ , where the angular 
brackets stand for averaging with respect to disorder 
configurations. The main order of {SQ'^) in l/kpl <C 1 
is given by perturbation theory diagrams shown in 
Fig[IIa))(b) [20|- These diagrams consist of Hikami boxes 
(shaded blocks in Fig. [Ha),(b)) connected by the wavy 
lines, which represent the sum of ladder diagrams: dif- 
fusons and Cooperons [l^]. Cooperons are strongly sup- 
pressed in magnetic fields in which magnetic fiux is larger 
then flux quantum per sample area. Since the UCF 
are usually studied experimentally in such a high mag- 
netic field regime, here we neglect the contribution of 
the Cooperon diagrams. In contrast, diffusons which are 
Green functions of quantum diffusion equation are not 
suppressed by a magnetic field although as we will show 
below their contribution depends on the efficiency of the 
symmetry-breaking disorder in the system. Below, we 
obtain diffusons by analyzing quantum kinetic equation 
in disordered graphene using the semiclassical approxi- 
mation (kpl ^ 1). 

Quantum kinetic equation describes relaxation of non- 
equilibrium inhomogeneous distribution of electrons in 
a disordered system. Using matrix Green functions 
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G ^ lVrG,TdtG. Advanced and retarded components 
of G are taken into account in the zeroth order of the 
gradient expansion, 



This determines the semiclassical Green functions 
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R/A 



±-(/ + Sn), 



where n = p/p- Dyson's equation for Keldysh com- 
ponent of disorder-averaged semiclassical Green function 
reads, 



FIG. 1: (a),(b) The diagrams which contribute to the main 
order in the diagrammatic expansion of the conductivity- 
conductivity correlation function. Here the solid lines stand 
for the impurity averaged retarded or advanced Green func- 
tions, the short wavy tails stand for the current vertices and 
the long wavy lines stand for the diffusion ladders. (c),(d) 
Hikami boxes of two types and additional diagrams which 
determine renormalization in the main order in l/kpl <^ 1. 
The dashed lines correspond to the disorder potential, (e) 
Diagrammatic equation for renormalized current vertex. 



in Keldysh representation, (^(ri, ii, ri/, ii/), we derive 
semiclassical form of Green functions for electrons in 
graphene. We separate slow and fast variables in G: 
r = ^(ri -I- ri') and t = ^{ti + ty), which vary 
at r 1/kp and t ^ fi/^p, and Sr — Vi — ri' and 
dt — (ti — ty), which vary at 6r ~ l/kp and St ~ h/ep. 
We then take the Fourier transform of 5G with respect 
to dt and Sr. The Fourier-transformed Green function, 



G(e,p,r, t) 



G«(e,p,r,i) G^(c,p,r,i) 
G^(e,p,r,t) 



obeys the following Dyson's equation. 
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S 



e-wSp-S" G = 1 



(yG(p,e,r,t)f), 



(2) 



where S is the self-energy matrix. 

The semiclassical Green function which describes low 
energy properties of the system is defined as 

g(e,n,r,t)= ^ j dCG{e,n{e + 0,r,t) 

where integration over ^ is performed in the vicinity of 
the Fermi level [2l|. The response of the system to ex- 
ternal perturbation is described by occupation numbers 
determined by the Keldysh component of semiclassical 
Green function. An equation for the latter can be de- 
rived using the gradient expansion, that is, assuming that 
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where 
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S-^ (/ + En) 0, 
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(3) 
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Note that, in Eq. ([3]) the energy e is defined with re- 
spect to mass surface shifted due to effects of disorder 
and Fermi line warping [23 |. 

Analyzing the main term in Eq. ([3]), eI]z(/+Sn)S25^, 
we find that the leading contribution to is propor- 
tional to the matrix (/ + Sn). Using 



5^'= E [9\i + ^rv)+Sg[ 



Ai 



(4) 



i—O.x.y.z 



where 5' are functions of n and I — Q,x,y, z, we consid- 
ered Sgl as a small correction (with an arbitrary matrix 
form 23]) and checked that the latter can be neglected 
in the leading order of the gradient expansion. 

Kinetic equation is obtained from equation ^ by sub- 
tracting its hermitian conjugate. After substituting the 
self-energies and g'' in the form ^ we find that. 



dtg' + «nv.g' + - (5' - (.9% - (5'n%n) 



-sHg'' 



-^'(.9'n%i 



0. (5) 



The angular brackets (...)<p in Eq. ^ and below denote 
averaging over momentum directions, and the coefficients 
T]' and (5' are defined as 



0, (5^ 
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The effect of the Fermi hne asymmetry is taken into ac- 
count by 



= (1-4^2), 



whereas Tw with (/, /') (x, y), {y, x) are equal to zero. 

The gradient expansion of Eq. ([5]) leads to the diffusion 
equation for the angle-average density matrix {g'')ip, 



where 24 



Di^v'4j2, r,V = 2r/(l + ryV), 



(6) 



Here, the valley-dependent transport times t^^ are de- 
termined by the efficiency of backscattering of electrons 
in corresponding mixed valley states (described by den- 
sity matrix components (i?'), {I = 0,x,y,z)). Such 
that, for example, due to the chirality of electrons in 
graphene backscattering off the potential disorder is sup- 
pressed [31, and in a sample with purely poten- 
tial disorder the transport time is given by r^^ — 2tq 
{I — 0,x,y,z). In realistic samples symmetry-breaking 
disorder restores the backscattering, which results in re- 
duced valley-state-dependent transport times r^^ < 2tq. 
The relaxation gaps F' are induced in this picture by 
symmetry-breaking disorder and the Fermi line warp- 
ing effect. This warping effect determines the different 
evolution operators of electrons in two different valleys, 
and suppresses the inter-valley coherence terms in the 
density matrix (g^/^). At the same time, however, the 
intra- valley components of density matrix (g^^^) are not 
affected by the trigonal warping. This effect is taken into 
account using a relaxation time, 

(we assumed t^^ ^ ''"^^)- Similarly, disorder terms 
mosAqEs and UzsAzS^, s — 0,x,y,z scatter electrons in 
different valley states differently, which also leads to re- 
laxation of (g^^^) without affecting (5"^^). Finally, inter- 
valley disorder terms Uxs^x'^s and Uysh^yYig mix the two 
valley states and lead to relaxation of all "valley-triplet" 
components of density matrix (g'^^*'^^), which is taken 
into account in Eq. ([6]) by the inter-valley relaxation rate 

Diffusons, can be now found as Green functions of 
diffusion equations ^ with initial inhoniogeneous dis- 
tribution 5q. To describe UCF in a small graphene 
sample we solve diffuson equations ^ with boundary 
conditions at current contacts = 0. The physical 
edge of graphene is atomically sharp and hence gen- 
erates strong inter-valley scattering, which suppresses 
valley-triplet diffuson modes near the edge thus leading 
to the boundary condition, D^/J'/^ = 0. In contrast, 
the particle-density ("singlet") mode has boundary 



condition (nV)!?" = corresponding to the absence of 
charge current through the edge. Solutions of the diffu- 
sion equations ^ for a rectangular graphene wire Lx 
are given by, 



I?'(r,r') 
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/II r \ (22./ riTTX \ I m-ny 
CmW = ^/— — sm — I cos f 

J^x J^y \ i-'x 



(7) 



where we take into account dephasing due to inelastic 
processes r^. 

As compared to the conventional electrons systems, 

Hikami boxes B[ — ^e^i^gt^r^ ^^[^ , B2 = je^vf^vr'^ ^^[^ 

and the current vertex Vx = vqt^j,/tYjx in monolayer 
graphene have to be renormalized by additional diagrams 
shown in Fig. [Hc),(d) and by vertex corrections (black 
dots in Fig. [1]). Both of these corrections contribute to 
the variance of conductivity in the main order in 1/kpl 
and are non-vanishing since current operator in mono- 
layer graphene is momentum independent 

The variance of conductance fluctuations is a sum of 
diagrams shown in Fig. [lja,b), in which the diagram 
Fig. [TJa) is encountered twice in the diagrammatic ex- 
pansion |25| and hence has a combinatorial pre- factor 2. 
As a result at T = we get jH], 
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It is interesting to compare UCF for rectangular phase- 
coherent graphene samples with ^ Ly and Lx ^ Ly. 
As mentioned before, we consider the system in an in- 
termediate magnetic field. In a narrow wire, V^^y^^ de- 
cay at the length ^ Ly, and the variance of conduc- 
tance is dominated by the "valley-singlet" diffuson com- 



ponent, 2?°, and, {60^ 



15 \ h J 



27| . which coin- 



cides with the standard result for quasi ID metallic wires 
in the unitary limit [20l |. In contrast, in the case of 
Ly 3> Lx, all diffuson components I?',^ = 0,x,y,z may 
contribute to the variance, depending on whether the 
effect of trigonal warping induces suppression of inter- 
valley diffuson components V^^^ or not. This determines 

{Sg^) = a^^t^ (t) (C(") is Riemann's zeta func- 
tion) with a — A for Lx < \J Dqt^ , a — Cq + Cz ^ 2 ioi 
^/DqtI < Lx < y/Dan, and a = 1 for y/DoT\i < Lx- 

Inelastic processes such as electron-electron or 
electron-phonon interactions limit the coherence length 
Lip ~ \/Dqt\p < Lx/y in the sample. In this case conduc- 
tance of the sample is determined as conductance of a 
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network of resistors of size each with conductance 
variance given by ([5]). Also, if at high temperatures 
Lt = KDq / kT < L and Lt < thermal broadening 
produces an additional self-averaging reducing the con- 
ductance fluctuations. For a square sample of size L x L: 

In conclusion, we have shown that the variance 
of interference-induced conductance fluctuations in 
graphene is of the order of the usual UCF value in metals, 
with a pre-factor dependent on the strength of the inter- 
valley scattering and shape of the sample. In a long wire 
of graphene or in the material with strong inter-valley 
scattering, the magneto-fluctuations of conductance have 
the variance typical for the unitary symmetry-class (in- 
termediate magnetic field). In a wide graphene sample 
(Lj. < Ly) with weak inter-valley scattering, the size of 
magneto-conductance fluctuations is increased as com- 
pared to unitary symmetry-class result by at least a fac- 



tor 2. This behavior is opposite to what was found 
for the weak localization magnetoresistance Q: the lat- 
ter was suppressed in the case of weak inter- valley scat- 
tering, whereas strong inter-valley scattering was found 
to restore the weak localization effect. This behavior 
contrasts the observation that in usual metals with non- 
chiral electrons UCF scale similarly to the weak local- 
ization correction to conductivity, made by Aleiner and 
Blanter [2^ . The analysis of the UCF in bilayer graphene 
showed a result very similar to the monolayer case despite 
a difference 0, in the electronic spectrum. 

The authors wish to thank E. McCann, B. Altshuler, 
A. Morpurgo, C.W.J. Beenakker, and A. Savchenko for 
helpful discussions. This project has been funded by EP- 
SRC grant EP/C511743 and ESF FoNE network SpiCo. 
We also thank K. Efetov for attracting our attention to 
that our UCF results for the intermediate asymptotic 
regimes coincide with those in [soj . 



[1] K.S. Novoselov et ai, Nature 438, 197 (2005); K.S. 
Novoselov et ai, Science 306, 666 (2004); Y. Zhang et 
ai, Phys. Rev. Lett. 94, 176803 (2005); Nature 438, 201 

(2005) . 

[2] K.S. Novoselov et ai, Nature Physics 2, 177 (2006). 

[3] A. H. Castro Neto et ai, cond-mat/0709.1163. 

[4] C.W.J. Beenakker, cond-mat/0710.3848. 

[5] E. McCann et ai, Phys. Rev. Lett. 97, 146805 (2006). 

[6] I.L. Aleiner and K.B. Efetov, Phys. Rev. Lett. 97, 236801 

(2006) . 

[7] K. Kechedzhi et ai, Phys. Rev. Lett. 98, 176806 (2007). 
[8] H. Suzuura, T. Ando, Phys. Rev. Lett. 89, 266603 (2002). 
[9] M. Titov, C.W.J. Beenakker, Phys. Rev. B 74, 041401(R) 

(2006) . 

[10] S.V. Morozov et ai, Phys. Rev. Lett. 97, 016801 (2006). 
[11] R.V. Gorbachev et ai, Phys. Rev. Lett. 98, 176805 

(2007) . 

[12] H. B. Heersche et al, Nature 446, 56 (2007). 

[13] F. V. Tikhonenko et ai, arXiv:0707.0140, 

[14] R. V. Gorbachev et ai, arXiv:0708.1700 

[15] A. Rycerz et al, Euro. Phys. Lett. 79, 57003 (2007). 

[16] Corners of the hexagonal Brillouin zone are K± — 

±(|/ia~^,0) (a is the lattice constant). 
[17] P.R. Wallace, Phys. Rev. 71, 622 - 634 (1947). 
[18] T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 



[19] Isospin and pseudospin matrices are defined as. 






1 










" 








1 













i 

















-1 








i 








-1 










-I 





1 













" 


1 








-1 








,Ax = 








-1 








1 





1 

















-1 







-1 











—i 


" 




1 

















i 







1 





i 














-1 








~i 
















-1 



This matrices determine two commuting subgroups of the 
group U4 of unitary transformations of a 4-component 
<I>: an isospin (sublattice) group SU2 = {e'""'^} and a 
pseudospin (valley) group SU^ = {e*'"'^}. In this basis 
time reversal T{W) of an operator is given by T{W) = 

AyEyWAyEy. 

[20] B. L. Altshuler JETP Lett. 41 648 (1985) [Pis'ma, Zh. 
Eksp. Teor. Fiz. 41, 530 (1985)]; P.A. Lee and A.D. Stone 
Phys. Rev. Lett. 55 1622 (1985); P. A. Lee, A.D. Stone, 
H. Fukuyama, Phys. Rev. B 35 1039 (1987). 

[21] J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 
(1986). 

[22] Integration of Eq. with respect to ^ gives rise to terms 
of the form ^ / d5^G(e,n(e -I- ^),r,t) = Aeg{e,ii,r,t), 
where Ae is a shift of the mass surface due to disorder and 
Fermi line warping. Such terms can be eliminated from 
the equation by recalculating the energy e with respect 
to the shifted mass surface. 

[23] Any matrix in the 4D isospin space can be ex- 
panded in the basis of linearly independent matri- 
ces [{i + En), (J + Sn)E,, E, (/ + En), E, (/ + Sn)E,]. 
Substitution of the solution in form of the Eq. Q into 
the main order in the Eq. ^ leaves only terms propor- 
tional to the last three matrices in the chosen basis in Sgl, . 



5 



[24] Conductivity of graphene monolayer with generic disor- 
der is given hy a — Ae^vDo = Ae^vv'^Ttr/'2 (where we 
took into account double spin degeneracy of carriers). 
Note that due to the valley-block-diagonal form of the 
current operator w^; = vT,^ only r^^ contributes to the 
conductance of graphene. At the same time the variance 
of conductance contains contribution from mixed valley 
states and therefore depends on r^V, / = 0, a;, y, z. 

[25] B. Altshuler, D. Khmelnitski, JETP Lett. 42, 359 (1985). 

[26] We take into account double spin degeneracy of carriers. 



[27] This result suggests that the anomalously large variance 
of sample-to-sample conductance fluctuations observed 
in ,153 is caused by a mechanism different from the phase 
interference which gives rise to the UCF. 

[28] 1. Aleiner and Ya. Blanter, Phys. Rev. B 65, 115317 
(2002). 

[29] E. McCann and V.I. Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006); J. Nilsson et ai, Phys. Rev. B 73, 214418 (2006). 
[30] M. Kharitonov and K. Efetov, arXiv:0801.0302. 



